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ABSTRACT 


This  paper  is  concerned  with  periodic,  laser- induced,  chemical  vapor 
deposition  recently  observed  experimentally.  In  order  to  inquire  further 
into  this  phenomena,  it  is  first  necessary  to  find  a  simple  means  of 
calculating  the  plasmon  field  strength  for  relatively  deep  gratings.  The 
Rayleigh  hypothesis  is  assumed,  and  only  p-po lari red,  normally  incident 
light  is  considered.  A  closed-form  equation  for  the  plasmon  field  intensity 
is  then  derived.  Also  discussed  is  the  behavior  of  the  plasmon  dispersion 
relation  for  a  shallow  grating,  but  for  a  complex  dielectric  constant  where 
the  imaginary  part  is  not  necessarily  small. 
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1 .  Introduction 
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Our  laboratory  has  been  investigating  laser- induced  chemical  vapor 
deposition  processes  as  recently  observed  experimentally.  In  particular. 

> 

|  wo  have  been  interested  in  the  periodic  structures  observed  by  Brueck  and 

i 

*  3 

Ehrlich  .  In  a  later  paper,  we  shall  discuss  this  experiment  directly. 

This  paper,  however,  will  be  devoted  to  developing  the  theoretical  framework 
necessary  for  the  calculations. 

The  periodic  structure  that  is  observed  evolves  from  the  plasmon  field 
induced  by  the  laser.  The  laser  interacts  with  the  random  surface 
roughness.  Because  the  plasmon  is  resonant  at  one  frequency  (for  shallow 
gratings),  that  frequency  component  grows  as 

|f  “  «I(x),  (1) 

where  £  is  the  grating  amplitude,  I(x)  is  the  total  field  intensity  as  a 
function  of  position,  and  a  is  the  proportionality  constant.  From  this  it 
is  clear  that  we  need  to  develop  expressions  for  the  plasmon  field  intensity 
as  a  function  of  grating  height  and  wavenumber.  We  can  simplify  the  problem 
considerably  in  all  that  follows  by  assuming  perpendicular  incidence,  in 
which  case  the  x-dependence  of  the  intensity  depends  only  on  the  plasmon 
field  and  not  on  the  reflected  or  incident  fields.  We  can  simplify  the 
problem  even  further  by  using  the  Rayleigh  expansion.  This  is  an  exact 
solution  above  the  selvedge  region  and  an  approximation  within  the  selvedge 
region.  However,  we  are  primarily  interested  in  field  strengths  above  the 
selvedge  region,  so  that  this  expansion  is  justified. 

Many  authors  have  discussed  the  interaction  of  radiation  with  a  rough 
surface  or  grating.  Petit^  has  written  a  good  introduction  to  the  subject. 
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w 

i 


Maradudin"  has  discussed  the  behavior  of  surface  polaritons  and  plasmons 

(which  together  are  what  we  are  generically  referring  to  as  a  surface 

plasmon).  In  particular,  he  has  worked  out  the  theory  for  randomly  rough 

surfaces  and  has  written  expressions  for  the  dispersion  relation  using  both 

the  Rayleigh  expansion  and  the  exact  formula  using  Green's  functions.  He 

shell  make  much  use  of  this  work  as  we  go  along.  Other  authors  have  also 

made  significant  contributions.  Jha,  Kirtley  and  Tsang^  have  developed  the 

theory  for  a  shallow  grating.  These  authors  have  made  extensive  use  of  the 

work  of  Toigo,  Marvin,  Celli  and  Hill,^'®  who  have  written  the  Rayleigh 

9  10 

expansion  in  a  remarkably  simple  form.  Finally,  Agassi  and  George  *  have 

developed  a  dressed  Rayleigh  expansion  which  eliminates  many  of  the 

numerical  difficulties  encountered  previously. 

We  are  interested  in  the  plasmon  field  intensity  above  the  selvedge 

region  for  reasonably  deep  gratings.  We  need  an  expression  as  a  function  of 

both  grating  height  and  grating  wavenumber.  Since  we  shall  also  need  the 

derivatives  of  this  function,  it  is  important  to  derive  an  expression  in 

closed  form.  To  date,  all  calculations  for  deep  gratings  have  involved  the 

inversion  of  large  matrices.  Given  the  conditions  of  our  application,  this 

is  impractical.  Therefore,  we  must  develop  an  approximation  which  is  still 

valid  in  the  deep  grating  case,  but  which  provides  a  closed  form  equation 

for  the  plasmon  field  strength.  Two  recent  papers  ’  have  made  progress 

12 

in  this  field.  In  particular,  the  work  by  Weber  represents  a 
generalization  of  our  present  result.  We  will  compare  our  result  with  his 
in  the  conclusion. 

We  should  also  point  out  that  there  has  been  considerable  work  done  on  a 

13 

closely  related  problem.  Fauchet  and  S legman  have  discovered  ripples  as  a 

result  of  laser  annealing,  and  Sipe  et  al.^  have  developed  an  elegant 


4 


theory  to  account  for  this  phenomena.  But  these  are  effects  of  surface 
damage  rather  than  surface  deposition,  and  hence  the  theoretical 
requirements  are  store  stringent  in  that  the  Rayleigh  hypothesis  can  no 
longer  be  applied. 

But  there  is  one  additional  problem.  The  experiment  of  Brueck  and 
Ehrlich  involved  a  plasmon  wave  along  a  cadmium  surface.  The  dielectric 
constant  for  cadmium  (at  UV  frequencies)  is  approximately 

e  -  -2.5  +  1.3i.  (2) 

Most  calculations,  even  for  shallow  gratings,  have  assumed  that  the 

4  5  6 

imaginary  part  of  the  dielectric  constant  is  either  zero  ’  .  or  small. 

Since  this  is  not  the  case  for  us,  we  shall  digress  a  bit  and  discuss  the 
effect  on  the  dispersion  relation  as  the  imaginary  part  of  the  dielectric 
constant  becomes  large. 

The  remainder  of  this  paper  is  organized  as  follows.  In  section  2  we 
shall  discuss  the  dispersion  relation  as  a  function  of  dielectric  constant. 
In  section  3  we  shall  derive  our  approximate  formula  for  the  plasmon  field 
strength.  Section  4  contains  a  discussion  of  our  results,  including  a 
numerical  comparison  of  our  approximation  with  the  exact  calculation. 

2.  Dispersion  Relation 

We  can  write  the  solution  of  the  homogenous  Helmholtz  equation  as 

m 

-a  "  EiPt.exP(itkx-o0z])  +  £  Atpl+exp(i[ktx+atzJ) 

t— • 

(3) 


*X(x) 


*<5(x). 


where 


k  *  5  c  ■  E1  ♦  lc2  h  <  •* 


kt  *  kt  ♦  tk»  «,  '  <‘2  -  kt2)1/2 


Ima  >  0 
°0  >  0 


(5) 


»,  -  <ek2  -  k,2)1/2 


Im6t  >0  Ptl  *  f[k(i  +  ati) 


and  where  k  is  the  plasnon  frequency  and  E.  the  incident  field  strength. 

o  * 

The  notation  has  been  culled  from  refs.  5  and  9.  This  solution  is  exact  as 

long  as  one  is  outside  the  selvedge  region.  Within  the  selvedge  region  it 

nay  be  a  good  approximation,  but  as  mentioned  previously,  we  are  prisiarily 

concerned  with  the  region  outside,  and  hence  we  can  use  eqs.  4  and  5  without 

reservation.  From  the  above,  it  is  apparent  that  Aq  is  the  reflected  field, 

while  A^  is  the  t-th  order  Bragg  reflection  (a^  real),  or  surface  plasmon 

(a,  imaginary).  We  are  interested  in  the  plasmon  effect,  and  hence  we  shall 

always  assume  that  k  >  k. 

8 

Now,  suppose  that  the  surface  profile  is  sinusoidal  and  can  be  written 

as 


5(x)  ■  £cos(k  x) 


(6) 


6 


and  that  £  is  small.  Then,  for  perpendicularly  incident,  p-polarized  light, 

6 


we  can  write  the  the  plasmon  field  strength  as 


Vl^at  -  fij] 


ek2  -  kg2(l+e) 


(7) 


Por  small  £,  is  very  small  unless  the  denominator  is  very  small.  If  e  is 
real,  the  denominator,  which  we  shall  call  the  resonance  factor  and  denote 
by  R,  is  zero  if 


2  2  1 

k  -  k 1  7-r  . 

g  e,  +  1 


(8) 


This  is  the  familiar  plasmon  resonance  condition  for  a  flat  surface  and  a 
real  dielectric  constant.  However,  if  t  is  complex  then  the  strength  of  the 
plasmon  field  depends  on  the  modulus  of  R,  and  it  is  no  longer  sufficient 
merely  to  minimize  the  real  part.  We  must  then  minimize 


| R | 2  -  IkJd+Cj)  -  k2El]2  +  te2(k2  -  k2)]2.  (9) 


Differentiating  with  respect  to  k  and  setting  the  result  to  zero,  we  obtain 

8 


a  new  dispersion  relation 


2  k^e^l+ej)  -  c2] 


g 


7  7 

(1+e^)*  +  e2 


(10) 


W 19  WWX"J»  PW*  l«n  W1V»  IffW  JTTTOimiWW.V, W, * 
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As  Cj  ■♦  0  this  reduces  to  eq.  S.  We  also  note  that  if  *  0,  then  |R|  *  0. 
and  the  resonance  is  both  dampened  and  broadened.  Equation  10  is  the  flat 
surface  plasmon  dispersion  relation  for  an  arbitrary,  complex-valued 
dielectric  constant. 

As  the  grating  becomes  deeper,  R  also  depends  on  higher  powers  of 
whereby  a  branching  occurs  in  the  dispersion  relation  and  the  resonance 
becomes  less  pronounced.  Thus  the  mechanism  of  the  periodic  deposition 
process  should  now  be  clear.  In  the  beginning,  the  incident  light  is 
resonant  with  one  frequency  component  of  the  randomly  rough,  but  very 
shallow  grating.  Because  of  this  resonance,  the  plasmon  of  that  frequency 
is  much  stronger  than  the  others.  Hence  the  field  intensity  is  periodic, 
and  deposition,  following  eq.  1,  has  the  same  periodicity,  leading  to  a 
growth  in  amplitude  of  the  grating. 

\ 

As  the  grating  becomes  larger,  the  dispersion  relation  branches,  and  the 
resonance  is  correspondingly  broadened  and  dampened.  This  results  in  a 
decrease  in  the  periodicity  of  the  field  intensity  and  a  corresponding 
decrease  in  the  differential  rate  of  deposition  (i.e.,  deposition  occurs 
evenly,  not  necessarily  more  slowly).  Hence  the  grating  stops  growing. 

3.  Calculation  of  Plasmon  Field  Intensity 

7  9 

We  can  write  an  expression  for  the  coefficients  of  eq.  3  as  * 


l  M..t\  - 


(ID 


where 


a  p,  +  k.k 
ml  t  m 


n-t, 


V,-  At  (1)  vt«K  •  i»» 

t  in 


(12) 


-an3  +  k-k 

U  "  . °.Tf!  <i>°  J<^°0  +  Pj>» 

m  aft  t  p  u  m 


(13) 


V. 

V, 

i 

>, 

> 

>v 

' 

>. 

I 


*0  "m 

and  J  (x)  is  the  n-th-order  Bessel  function.  In  general,  this  results  in  an 

n 

infinite  number  of  coupled  equations,  from  which  the  A^'s  can  be  extracted. 

9 

In  practice,  Agassi  and  George  have  shown  that  50  equations  is  sufficient 
to  find  an  accurate  solution.  However,  in  our  circumstance  this  is 
impossible.  We  must  find  the  plasmon  field  strength  as  a  function  of  both 
grating  height  and  wavenumber,  which  means  that  we  would  have  to  invert  two 
50  x  50  real  matrices  for  each  data  point.  Using  the  approximation  scheme 
which  we  shall  derive,  we  have  calculated  10,000  data  points.  This  is 
perhaps  more  than  necessary,  but  doing  so  has  enabled  us  to  thoroughly 
understand  the  behavior  of  these  functions. 

For  our  approximation,  we  borrow  a  page  from  Maradudin^  and  extend  that 
result  to  higher  orders.  The  method  rests  on  two  ideas.  One  is  that  if  A^ 
is  a  plasmon  at  resonance  with  k^,  then 


A.  »  A 
1  m 


m  *  1. 


(1A) 


K 


K 


The  second  condition  is  to  note  that 


M  -  1,  M  M  -  C, 

nn  n,n±l  n,n±2 


(15) 


so  that  even  for  reasonably  deep  gratings  we  can  assume  that  > 


etc. 


A  p^ecial  case  of  eq.  11  is 


(16) 


MnAi  + 


)  M,  A  -  m,E4 

i  1,0  m  H1  i 


o*l 


or 


M,  E, 


X* 


m*l  In  n 
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(17) 


Thus  R  *  M^,  which  oust  therefore  go  to  zero  if  is  to  be  resonant. 
Expanding  A^  to  first  order  in  £  yields  our  result  of  eq.  4. 

Equation  17  is  exact.  Now  we  nake  our  approx imat ion.  Suppose  that  n  * 
1.  Then  we  can  write 


M  A 
no  m 


l 


M  A 

mp  p 


*oEi 


p*n 


(18) 


But  because  we  are  at  resonance,  of  the  sun  over  p,  p*l  is  by  far  the 
largest  tern.  Thus  we  can  sinplify  eq.  18  as 


V*i 


,n 


M 


(19) 


Substituting  eq.  19  into  17,  we  obtain 


In  the  case  of  normal  incidence,  ■  A  ^  by  symmetry,  and  we  are  not 
Justified  in  leaving  A out  of  eq.  19.  Inclusion  of  this  term  for  normal 
incidence  yields 


A 


1 


<ui  - 1  <HiP  *  "iVs2  )1!i 
.  in  n _ 

[<MU  7  M,.,)  -  l  (Hlp+  M^HN^ 

P*1  PP 


(21) 


where  p  is  summed  only  over  non-negative  integers.  It  is  clear  that  in  the 

shallow  grating  case,  both  eqs.  20  and  21  reduce  to  eq.  7. 

We  shall  use  eq.  21  with  one  modification:  To  ensure  numerical 

9 

convergence,  we  shall  use  the  dressed  Rayleigh  expansion.  Thus  we  rewrite 
eq.  10  as 


(22) 


where 


m,n 


M  expUCU  +  2. 

O)  n  n  n 


i) 


(23) 


uexpdB  O  , 


A^  ■  A  exp(ia  O* 


This  eliminates  numerical  instabilities,  not  only  in  the  exact  calculation 
but  also  for  our  approximation.  In  practice,  for  the  grating  heights  of 
Interest,  we  have  found  that  summing  p  from  0  to  15  in  eq.  21  provides 
sufficient  accuracy. 


4.  Conclusion 


Equation  21,  modified  if  necessary  by  the  transformation  of  eq.  23, 
constitutes  our  closed  form  approximation  for  the  plasmon  field  strength. 

Let  us  consider  the  advantages  and  shortcomings  of  our  approximation.  The 
first  constraint  is  to  notice  that  for  our  approximation  to  be  valid,  eq.  14 
must  hold.  This  means  that  we  must  be  on  the  resonance  frequency.  Were  we 
not  on  the  resonance  frequency,  then  the  approximation  used  in  deriving  eq. 
19  would  not  be  valid.  Secondly,  our  approximation  still  is  a  function  of 
grating  height.  Note  that  we  have  used  eq.  15,  which  means  that  we  are 
discarding  higher  orders  of  £  for  non-resonant  plasmons.  As  these  become 
large,  then  our  approximation  becomes  less  accurate. 

Nevertheless,  we  are  happy  to  report  that  this  method  preserves  the 
qualitative  features  of  the  full  calculation.  Figure  1  illustrates  this 
very  clearly.  Here  we  have  calculated  the  plasmon  field  strength  (|  Aj)  for 
various  grating  heights  using  both  our  new  method  and  the  full  calculation. 
It  is  clear  that  for  gratings  up  to  £  *  25  nm,  the  simple  calculation 
preserves  the  qualitative  features  of  the  exact  method.  This  corresponds  to 
a  value  of  k^£  a  .7,  which  is  a  very  deep  grating.  The  quantitative  picture 
stay  also  be  better  than  appears  in  the  graph.  The  calculation  was  done  for 
kg  ■  2.95  x  107  m  *.  This  is  the  resonance  frequency  for  shallow  gratings, 
but  it  is  off  resonance  for  deeper  gratings.  At  25  nm,  we  have  observed 
that  the  resonance  frequency  is  closer  to  3.1  x  10  7  m  *.  We  have  ignored 
this  shift  in  frequency  in  our  calculations. 

As  mentioned  previously,  our  result  compares  with  the  more  general 
12 

result  of  Weber.  His  is  more  general  in  that  it  explicitly  accounts  for 
three  terms  in  the  Rayleigh  expansion  rather  than  just  one.  However,  at 
resonance,  the  other  two  terms  can  be  expected  to  contribute  negligibly,  and 


hence  our  result  should  b«  nearly  as  accurate.  Off  resonance,  of  course, 
the  other  terms  are  significant  and  Weber's  method  is  clearly 
preferable.  Strict  numerical  comparison  is  impossible  since  he  did  his 
calculation  using  silver,  with  a  small  C2*  whereas  ours  is  for  cadmium. 
Furthensore,  his  example  is  for  a  sawtooth  grating  while  ours  is  for  a 
sinusoidal  grating.  But  it  is  very  clear  from  a  comparision  of  figure  1 
with  his  table  I  that  the  same  qualitative  features  hold,  nasMly  that  the 
plasmon  resonance  intensity  is  underestimated.  Our  formalism,  while  not  as 
general,  is  derived  much  more  simply  and  is  easier  to  use,  and  is  probably 
just  as  accurate  where  applicable. 

In  a  later  paper^  we  shall  apply  this  formalism  to  the  problem  of 
periodic,  laser- induced,  chemical  vapor  deposition.  We  should  be  able  to 
predict  the  growth  rate  and  maximum  peak  height  using  the  above  method. 
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Figure  1. 


FIGURE  CAPTION 


This  graph  compares  the  plasmon  field  strength  ( j | ) 
calculated  from  eq.  21  (using  the  dressed  Rayleigh  expansion) 
with  that  from  the  exact  method.  The  dielectric  constant 
used  is  -2.5  +  1.31,  the  incident  light  has  a  wavelength  of 
257  nm  (k  »  2.44  x  107  a  ^ ) , and  the  grating  wavenumber  is 
taken  as  2.95  x  107  m”*. 
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